Testing intermediate- age stellar evolution models with VLT 
photometry of LMC clusters. I. The data 1 



Carme Gallart 2,3 , Manuela Zoccali 5 , Giampaolo Bertelli 6,9 , Cesare Chiosi 7 , Pierre 
Demarque 4 , Leo Girardi 8 , Emma Nasi 9 , Jong-Hak Woo 4 , Sukyoung Yi 10 

CN " 

O ■ ABSTRACT 
O 
(N 

This is the first of a series of three papers devoted to the calibration of 
^ . a few parameters of crucial importance in the modeling of the evolution of 

intermediate-mass stars, with special attention to the amount of convective core 
overshoot. To this end we acquired deep V and R photometry for three globular 
clusters of the Large Magellanic Cloud (LMC), namely NGC 2173, SL 556 and 
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NGC 2155, in the age interval 1-3 Gyr. In this first paper, we describe the 
aim of the project, the VLT observations and data reduction, and we make 
preliminary comparisons of the color-magnitude diagrams with both Padova and 
! Yonsei-Yale isochrones. Two following papers in this series present the results 
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of a detailed analysis of these data, independently carried out by members of 
the Yale and Padova stellar evolution groups. This allows us to compare both 
sets of models and discuss their main differences, as well as the systematic 
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effects that they would have to the determination of the ages and metallicities 
of intermediate-age single stellar populations. 

Subject headings: stars: evolution, color-magnitude diagrams, galaxies: 
individual (LMC), clusters: individual (NGC 2173, SL 556, NGC 2155) 

1. INTRODUCTION 

In an epoch of extraordinary discoveries on the high-z Universe, we still have gaps in 
our understanding of stellar evolution, which is in fact essential for the correct interpretation 
of the light of distant galaxies. Galactic star clusters have traditionally provided a major 
way to study stellar evolution. However, our Galaxy only contains stars in age and chemical 
composition domains which reflect its particular history, and these are the domains which 
have been relatively well explored. Extrapolation outside the limits of the well explored age 
and metallicity range is not always safe. For instance, it has become increasingly apparent 
in recent years that metal content can affect stellar evolution in unexpected ways. The UV 
upturn in elliptical galaxies is a case in point. Studies of very metal rich stellar systems 
revealed that, contrary to simple expectations, old metal rich stellar populations do not 
simply become redder in all wavebands as they evolve, but rather produce a population 
of UV bright stars (Greggio & Renzini, 1990; Horch, Demarque & Pinsonneault, 1992; 
Fagotto et al. 1994; Yi, Demarque & Oemler 1998). On the other hand, understanding the 
evolution of extremely metal-poor and metal-free (Population III) stars will similarly be 
essential for the interpretation of primordial, high-z stellar populations. 

In intermediate-age stellar populations, the color-magnitude diagram (CMD) and the 
luminosity function are affected by convective core overshoot. While there is not a general 
agreement on the efficiency of this process, and its extent is not well established, most 
researchers agree that this parameter affects significantly the morphology of the CMD 
of these clusters and naturally the determination of their ages as well (e.g. Rosvick & 
VandenBerg 1998; Keller, Da Costa & Bessell 2001; Meynet, Mermilliod & Maeder 1993; 
Carraro et al. 1993; Demarque, Sarajedini & Guo 1994). This uncertainty is thus also a 
problem for the spectral dating of distant stellar systems from their integrated light (Heap 
et al. 1998; Yi et al. 2000). In addition, deep CMDs of intermediate-age clusters are 
essential to help disentangle the relative importance of other poorly determined parameters, 
such as mass loss during red giant branch (RGB) evolution, internal rotation (important 
for the more massive objects), and helium content. In particular, the adopted value of the 
parameter 5Y/ 5Z of helium enrichment significantly affects the mass luminosity relation of 
stellar models, and consequently their evolutionary lifetimes. 
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The papers in this series deal mainly with the problem of convective core overshoot in 
intermediate-age stellar populations, and therefore we will briefly review some of the recent 
work in the subject. Early observational arguments in favor of convective core overshoot 
come from clusters like Pleiades (Maeder k Mermilliod 1981; Mazzei k Pigatto 1989) and 
from other Galactic clusters, as discussed in great detail by Maeder k Mermilliod (1981) 
and Mermilliod k Maeder (1986). Barbara k Pigatto (1984) and Chiosi k Pigatto (1986) 
argued for overshoot in stars with masses in the range 1.5 — 2.2M by pointing out that 
the base of the RGB is not well populated in clusters with age 1-2 Gyr (whereas it is 
in older clusters), as if degenerate He- ignition and He-flash were avoided for this mass 
range, in contrast with classical models. More recent studies include those of Aparicio et 
al. (1990), Carraro et al. (1993), Meynet et al. (1993); Rosvick k VandenBerg (1998); 
Demarque et al. (1994); Dinescu et al. (1995); Kozhurina-Platais et al. (1997) which deal 
with Galactic open clusters of ages 1.0 to 6.0 Gyr. All conclude that a certain amount of 
convective core overshoot is preferred to reproduce the CMDs. Being Galactic, all these 
clusters have metallicity close to solar, and therefore do not allow the possibility of testing 
the dependence of convective core overshoot on metallicity. 

The Magellanic Clouds (MC) offer an unusual opportunity to test stellar evolution 
since populous clusters of a wide variety of ages and chemical compositions can be observed 
and compared at effectively the same distance. In particular, they contain young and 
intermediate-age metal-poor stellar populations that are absent in the Galaxy. A number of 
young MC clusters (age < 500 Myr) have been used by Lattanzio et al (1991), Vallenari et 
al. (1991, 1994), Stothers k Chin (1992), Brocato, Castellani k Piersimoni (1994), Chiosi 
et al. (1995), Testa et al. (1999) and Barmina et al. (2002) to constrain evolutionary 
models of young stars. It is interesting to note that these different studies do not always 
agree on the preferred amount of convective core overshoot to best reproduce the observed 
CMDs. In particular there is a long lasting debate about the young LMC cluster NGC 
1866, where turnoff stars have about 4 — 5M Q : Chiosi et al. (1989a,b) and Brocato k 
Castellani (1988), Lattanzio et al. (1991) and Brocato et al. (1994), Testa et al. (1999) and 
Barmina et al. (2002), are pairs of papers presenting systematically opposite conclusions 
(for or against overshooting). To briefly summarize only the most recent ones, Testa et al. 
(1999) obtained deep photometry (V < 24) of a wide region of this cluster, and concluded 
that the best fit to both the integrated luminosity function and the magnitude of the 
horizontal branch is achieved by classical models, without overshooting, provided that a 
30% fraction of binaries is allowed. However, a new analysis of the same data with a more 
accurate completeness correction and normalization, by Barmina et al. (2002), showed that 
models with overshoot give a better fit to both the overall morphology of the CMD and to 
the integrated luminosity function of main sequence stars, reproducing the correct ratio of 
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main sequence to post-main sequence stars. 

Thanks to the brightness of the stars in the young MC clusters, the above mentioned 
observations were feasible using medium size telescopes. Instead, our aim is to explore 
the less studied regime of intermediate-age, low metallicity clusters, which turnoff is 
by definition fainter, and therefore only 8-meter size telescopes under excellent seeing 
conditions, or the HST, can give the necessary high quality photometry. Our intention 
is also to independently analyze the data with two of the most recent sets of stellar 
evolutionary models widely in use, namely, those of Padova (Girardi et al. 2000) and 
Yonsei-Yale (Yi et al. 2001, Y 2 thereafter). 

In this paper, we discuss the project setup (Sec. 2), and we present the data obtained 
with VLT: the observations, photometry and crowding tests are presented in Sec. 3, while 
the procedure used to statistically substract the LMC field stars from the cluster CMDs is 
discussed in Sec. 4. In Sec. 5, we briefly discuss the main differences between the stellar 
evolution models of Girardi et al. (2000) and Yi et al. (2001). In Sec. 6, a preliminary 
comparison of the cluster CMDs with these stellar evolutionary models is performed. In two 
forthcoming papers that follow in this same issue, these data are independently analyzed in 
detail with both the Y 2 (Woo et al. 2002) and Padova (Bertelli et al. 2002) models, in such 
a way that some feedback on the stellar evolution models used can be provided. 



2. OBSERVATIONAL GOAL AND TARGET SELECTION 

The goal of this project was to obtain accurate V, R photometry down to about 
4 magnitudes below the main sequence turnoff of a number of intermediate-age MC 
clusters, in order to construct CMDs and luminosity functions to be used as test for stellar 
evolutionary models computed independently by members of the Yale Group (Demarque, 
Yi, Woo)and the Padova Group (Bertelli, Chiosi, Girardi, Nasi). 

We originally proposed to observe eight clusters, in pairs of 1, 2, 3 and 5 Gyr each 
(according to the ages quoted in the literature) to be able to check the extent of the 
fluctuations in the properties of clusters (e.g. luminosity function) of the same age. The 
candidates were selected according to: a) their age (from Geisler et al. 1997; Sarajedini 
1998; Mighell, Sarajedini & French 1998), b) their position in the galaxy, trying to avoid 
areas densely populated by field stars, and c) their richness. The original list of clusters 
included NGC 2209 and NGC 2249 (~ 1 Gyr old), NGC 1651 and NGC 2162 (~ 2 Gyr 
old), NGC 2173 and SL 556 (~ 3 Gyr old) and NGC 2155 and Kron 3 (~ 5 Gyr old). The 
first four younger clusters were to be observed with NTT down to a limiting magnitude of 
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V=23.5 and 24.5, while the rest were proposed for F0RS1@VLT, with the goal of reaching 
limiting magnitudes V=25.5 (for NGC 2173 and SL 556), V=26.0 for NGC 2155 and 
V=26.3 for Kron 3. Due to the high resolution needed to obtain accurate photometry in 
such rich and distant clusters, we requested service mode observations imposing a seeing 
constraint of FWHM <0.6". In the end, due to observational problems and restrictions, no 
observations with NTT were performed for this program, while 3 out of the 4 clusters to be 
observed with VLT were indeed observed, namely, NGC 2155, NGC 2173 and SL 556. 

In what follows we summarize the results of previous studies of these clusters, 
concentrating in particular on the metallicity determinations, being this an important 
input parameter for our analysis. The distance modulus of the LMC, being the first step 
of the distance scale ladders, has been the subject of an intense debate in recent years and 
trying to summarize it goes beyond the scope of this paper: the generally accepted value 
of (m — M) v = 18.5 (Benedict et al. 2002) has been adopted, allowing a variation of ±0.2 
in the isochrone fit, which can be justified by volume effect of the LMC on the distance 
modulus of each cluster. Small variation of the reddening around E(B — V) ~ 0.04 along 
the different lines of sight of the three clusters, will be allowed in our analysis as well. 

Available determinations of the metallicity of these clusters are very inhomogeneous, 
first of all because of the different methods applied, and second because some of them 
lead to the determination of the global metallicity, [M/H]—logZ /Z Q — logX/X Q (e.g. 
when using isochrone fitting), while some others give a measure of the iron abundance, 
[Fe/H]=logFe/Fe® — logX/X & (e.g. spectroscopic or photometric methods calibrated 
against [Fe/H]). Obviously the two quantities are different if the alpha-element ratio of 
the cluster stars is different from solar, as it is for galactic globular clusters. Since there 
currently is no high dispersion spectral analysis of the MC clusters from which one could 
infer their alpha-element enhancement, in what follows we will keep the distinction among 
the two chemical abundances, as measured by the various authors, reminding the reader 
that if the alpha-element enhancement of these LMC clusters were the same as that of 
galactic globular clusters (i.e., [a/Fe]~ 0.3, Carney 1996) then their global metal abundance 
[M/H] can be assumed to be ~ 0.21 dex higher than the iron content [Fe/H] (Salaris, Chieffi 
& Straniero 1993). 

Low resolution spectra of four stars in NGC 2173 were obtained by Cohen (1982), 
who found a global metallicity of [M/H]=-0.75. A lower value ([Fe/H]=-1.4) was obtained 
by Bica, Dottori and Pastoriza (1986) with integrated photometry, while the Ca triplet 
analysis of two stars, by Olszewski et al. (1991), gave [Fe/H]=-0.24. B, V photometry 
for this cluster has been published by Mould, Da Costa & Wieland (1986), who estimate 
[M/H]=-0.8 from the RGB slope. They found an age of either 1.4 or 2 Gyr, when adopting a 
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long (18.7) or short (18.2) absolute distance modulus, respectively. More recent photometry 
for NGC 2173 has been published by Corsi et al. (1994), who were looking for global 
properties of intermediate-age populations and therefore do not analyze the CMD of this 
cluster individually, but only superimposed to that of 10 other clusters. Finally, Geisler 
et al. (1997), adopting the metallicity [Fe/H]=-0.24, as found by Olszewski et al. (1991), 
estimate a cluster age of 2 Gyr, in agreement with that found by Mould et al. (1986). 

SL 556 (Hodge 4) is a relatively less studied cluster. No spectroscopic studies have 
been carried on so far. U, B, V photometry has been published by Mateo & Hodge (1986) 
who estimate a metallicity of [Fe/H]=-0.7 using various photometric indexes. Olzewski 
et al. (1991) determined [Fe/H]=-0.15, while more recent photometric studies have been 
published by Sarajedini (1998) and Rich et al. (2001), who analyzed the same HST data, 
and, adopting a metallicity [M/H]=-0.68 both found an age between 2 and 2.5 Gyr. 

NGC 2155 has been studied by Bica et al. (1986) who performed DDO integrated 
photometry of a number of MC clusters. They estimated a metallicity of [Fe/H]=-1.2. More 
recently, a significantly higher value ([Fe/H]=-0.55) has been found by Olszewski et al. 
(1991) who performed Ca II triplet spectroscopy of three cluster giants. Although possibly 
more accurate than the previous integrated photometry, this method is not very accurate 
because its calibration in the metal-rich ([Fe/H]> —1), relatively young populations is more 
uncertain than for the old, metal poor ones (Pont et al. 2002). Indeed, for the three clusters 
studied here, the metallicities quoted in Olszewski et al. (1991) are systematically higher 
(by ~ 0.5) than any other published value. A more recent determination of the cluster 
metallicity from the RGB slope comes from the HST photometry published by Sarajedini 
(1998). The latter estimates [Fe/H]=-l.l, and an age of ~ 4 Gyr. The same HST data 
used by Sarajedini have been re-analyzed by Rich et al. (2001) and compared with the 
isochrones by Girardi et al. (2000), giving a global metallicity of [M/H]=-0.68 and an age 
of 3.2 Gyr. 

3. THE DATA 

The three LMC clusters NGC 2173, SL 556 and NGC 2155 were observed with 
FORSl@VLT-UTl, in service mode, through the filters V and R. Figure [TJ shows a median 
of all the frames for SL 556, as an example of the image quality of the data. FORS1 was 
used in normal resolution: its pixel size of 0.2 "/px allows a good PSF sampling, while the 
relatively large field of view (6.8 x 6.8 arcmin) allowed us to sample, in the field outermost 
regions, a good portion of the background LMC field. Total exposure times and average 
seeing conditions for each cluster are listed in Table 1. 
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De-biasing and flat-fielding were performed with standard IRAF packages, using 
the calibration frames (bias and sky flat fields) obtained during the same night of the 
target observations, or one night before. Photometric reduction was carried out with the 
DAOPHOT II/ALLFRAME package (Stetson 1987, 1994). Stellar PSF models for each 
frame were obtained by means of a large set of uncrowded and unsaturated stars. After 
performing aperture photometry on each frame, spatial transformations were determined 
among all the frames of each cluster, in order to be able to register them and obtain a 
median image free from cosmic rays and having high S/N ratio. This image was used 
for star finding, in order to obtain a deeper catalog to use as input for the simultaneous 
profile-fitting photometry of all the frames, made with ALLFRAME. 

According to the log file of the service observations, the present data were obtained 
under unstable photometric conditions: only part of the nights were photometric, and 
some cirrus seemed to be always present close to the horizon. Therefore, calibration to the 
Johnson standard photometric system was performed by comparison of the stars in common 
with another set of images of the same clusters, obtained by our group at the CTIO 4m 
telescope using the MOSAIC camera, on the nights Nov. 30th to Dec 1st 1999, which 
were devoted to a different project (Gallart et al. 2002). Those images, although not deep 
enough for our purpose, were obtained under photometric conditions, and were already 
reduced and calibrated, by means of a set of Landolt (1992) standard fields observed on the 
same nights. The V, R cluster images were obtained on chip#2 of the MOSAIC camera, on 
which we observed around 50 Landolt standards on the color range (V — R) — 1. No 
color terms were considered necessary on the photometric transformation, and zero-points 
had typical errors of ± 0.004. Aperture corrections for the cluster stars were obtained from 
a list of ~ 70 relatively isolated, bright stars in each cluster, with standard errors of ±0.002 
mag. Therefore the total maximum zero-point in the calibration of the CTIO images is of ± 
0.005 magnitudes. The FORS1 data were then calibrated by means of the ~ 100 brightest 
and most isolated stars in common with the CTIO data, from which we derived a zero point 
and a color term, with a total error (including the one quoted above) of ±0.006. 

In the two following papers in this series, we will use both the morphology and the 
number of observed stars in different regions of the CMDs as compared with theoretical 
predictions. To this end, it is crucial to correct the observed counts for completeness. 
Completeness corrections were determined by means of artificial-star experiments. A total of 
about 30,000 stars were added to the original frames, with magnitudes and colors consistent 
with the main branches of the cluster and field CMD. In order to avoid overcrowding, while 
optimizing the CPU time, the artificial stars were added along the corners of an hexagonal 
grid, as explained in Zoccali et al. (2000), in five independent experiments. As usual, 
photometry of the artificial frames was performed in a way identical to that of the original 
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ones. 



4. UNVEILING THE CLUSTER CMDs: SUBTRACTION OF THE LMC 

FIELD STARS 

Figure | shows the CMDs of the three LMC clusters NGC 2173, SL 556 and NGC 2155. 
All the stars measured with good precision in the whole FORS1 frames are plotted here: 
therefore the CMDs include both cluster and field stars. In order to obtain a cleaner cluster 
CMD, we proceeded in the following way. 

For each cluster, an annulus around the center was selected with the criterion of being 
large enough to contain most cluster stars, but small enough to minimize the background 
contribution. The CMD for the stars in this annulus are shown in the upper left panel 
(panel a) of Figs. |3| to [| together with the inner and outer radii of the selected annuli. 
The small circular region containing the very center of the cluster was excluded because 
stars in these areas were too crowded, and therefore their photometry was rather poor. A 
CMD typical of LMC field stars was extracted from a second region, outside a given radius 
around the cluster center, selected with the criterion of containing a number of stars large 
enough to insure good statistics, but far enough from the cluster not to contain a significant 
fraction of cluster stars. The CMDs extracted from these "field" regions are shown in the 
upper right panels (b) of Figs. to [5[ together with the inner boundary of the region. The 
outer boundary is set by the edges of the frame. 

For each star in the field CMD (panel b) we picked up the closest (see below) star in 
the cluster CMD (panel a), and decided whether to subtract it or not according to a factor 
given by the ratio of the areas of the two regions, multiplied by the ratio of the completeness 
factors for the two regions (being less crowded, the field CMD is more complete than the 
cluster one). The distance, d, from the field star to each cluster star on the CMD was 
actually computed by enhancing the difference in color by a factor of 7 with respect to the 
difference in magnitude, i.e.: 

d= 07 x A(V - R)] 2 + AV 2 . 

which defines an ellipse with major axis along the magnitudes and minor axis along the 
colors. The reason for this is that given a certain star in the field CMD, its "twin" in 
the cluster CMD has a color that is better constrained than the magnitude, because any 
small difference in either distance, reddening, or mass between the two would result in 
a difference in magnitude larger than that in color. If there were no stars closer than a 
maximum distance of d = 0.42 magnitudes then no star was subtracted from the cluster 
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CMD, and the field star was flagged as unsubtracted. Both the enhancement factor and the 
maximum distance are somewhat arbitrary, and they were chosen after various experiments, 
according to two criteria: i. The CMD of the subtracted stars (panel d) should be as similar 
as possible to that of the field stars (panel b) and ii. the number of field stars flagged as 
unsubtracted should be as small as possible. 

The result of this procedure (i.e., the CMD in panel c) is indeed satisfactory, since in 
general only the CMD of a single intermediate- age population is present in this plot. Some 
younger main sequence stars are still present in the decontaminated CMD of NGC 2155, 
due to the fact that, being the young LMC population intrinsically clumpy (like any young 
stellar population in galaxies), the assumption we make here, that the stars we measure 
just outside the cluster boundary are representative of a population identical to that 
contaminating the cluster CMD, is true only as a first approximation. 

5. A COMPARISON BETWEEN THE YONSEI-YALE AND PADOVA 

STELLAR EVOLUTIONARY MODELS 

The goal of this project is to compare the observed and theoretical CMD of 
intermediate- age, metal-poor clusters, in order to constrain some of the still uncertain 
parameters used in stellar evolutionary models. In the following two papers of this series, 
the data presented here will be independently analyzed by members of the Yale and the 
Padova stellar evolution groups, with the aim of calibrating both sets of models, and also 
to better identify the effects of different assumptions on the simulated CMDs. In view of 
this, here we present a comparison between the sets of isochrones published by the two 
groups, in the range of ages and metallicity spanned by the observed clusters, namely l-j-3 
Gyr and Z=0.004-=-0.008. We consider this as a "first-order" comparison because we are 
considering only the published models, without any fine tuning of the input parameters, 
that will be performed in the following two papers. Because the two grids of models and 
isochrones were derived independently, a comparison between them can provide a measure 
of the uncertainties in isochrone fitting, as well as the robustness of the ages derived for 
LMC star clusters. 

Figure |B| shows a set of isochrones of ages 1, 2 and 3 Gyr, and Z=0.004, obtained from 
Girardi et al. (2000) and Yi et al. (2001) in the theoretical [log L/L Q ,log T e fj] and in the 
observational [V, (V — R)] plane. The two sets of isochrones are quite close to one another 
up to the main sequence turnoffs in the theoretical plane. Above that, the Y 2 isochrones are 
slightly brighter than the Padova ones, while the RGB is slightly cooler. Another apparent 
difference is in the 3 Gyr old isochrone, which is calculated from models with a moderate 



-10- 



amount of core overshooting (gradually increasing with mass in the range between 1.0 
and 1.5M , see Girardi et al. 2000) in the Padova set, but with no overshooting for Y 2 
isochrones by Yi et al. (2001), who adopted some overshooting only for isochrones younger 
than 3 Gyrs. The differences are more important in the observational plane (right panel), 
where the uncertainties in the color transformations add up. In the following we will discuss 
in some detail the differences in the input physics and color transformations between both 
sets of models, which may explain the differences apparent in the figures. 

In the theoretical log L/L Q vs. log T e ff plane shown in the left-hand panel of Figure H 
the differences must be due to differences in the stellar models used in constructing the 
isochrones. Regarding the input physics, both sets of isochrones make use of the OPAL 
opacities (Rogers & Iglesias 1995) for the stellar interiors, and of the Alexander & Ferguson 
(1994) opacities at low temperatures. But they differ in the choice of the equation of state 
in the interior. While the Y 2 isochrones are based on the OPAL equation of state (Rogers, 
Swenson & Iglesias 1996), the Padova isochrones use the prescription by Straniero (1988) 
and the MHD equation of state (Mihalas, Dappen & Hummer 1988). Also different are the 
energy generation rates and the details of the solar calibration, which determine the choice 
of two important input parameters, the initial helium abundance and the mixing length in 
the convective zone. The mass-luminosity relation for stars depends significantly on the 
chosen helium abundance, and the mixing length parameter in the convection zone of cool 
stars affects the radius of cool star models. 

Within the conventional method of overshoot parameterization Padova isochrones 
differ from Y 2 ones in the adopted constraints to the overshoot efficiency, as described 
by Bressan et al. (1981) and Girardi et al. (2000), and by Yi et al. (2001) and Woo & 
Demarque (2001). In the Y 2 models, the extension of the convective motions beyond the 
formal boundary of the convective core defined by the Schwarzschild criterion is measured 
from the convective core boundary outward (i.e. above the convective boundary). Its 
magnitude is expressed as a fraction of the local pressure scale height H p at the convection 
boundary. In the Y 2 models, OS=0.2 (meaning 0.2H p ) is adopted for younger isochrones 
(< 2 Gyr) and OS=0.0 for older ones (> 3 Gyr). In the Padova models, the overshoot 
length A c is measured across the formal Schwarzschild convective boundary (i.e. it is 
measured from inside the convective core and it straddles the convective boundary). The 
overshooting parameter A c is so defined that overshoot of A c = 0.5 (Padova) corresponds 
to approximately OS = 0.25 H p in the Yale code. Understanding this difference will be 
important to correctly interpret the conclusions on the need of a particular amount of 
overshoot in the models, in the next two papers in this series. 

Despite these main differences and some minor others in the construction of the stellar 
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models, it is encouraging to note the generally good agreement in the log L/L Q vs. log T e ff 
plane, and in the derived ages. 

In the CMD, shown on the right-hand side in Figure |^, the differences between the Y 2 
and Padova isochrones are more pronounced. This emphasizes the still large uncertainties 
in the color transformation, used to translate the isochrones from the theoretical log L/L & 
vs. log T e ff plane into the observational My vs. (V — R) plane. Bertelli et al. (1994) and 
Girardi et al. (2000) isochrones are based on the Kurucz (1992) library of stellar spectral 
energy distributions, convolved with filter spectral responses (extended at high temperature 
with black-body spectra and at low temperature with empirical M stars spectra as described 
in those paper). 

The Y 2 isochrones used here are instead based on the semi-empirical Lejeune, Cuisinier 
& Buser (1998) color transformations (BaSeL colors). The spectral libraries assembled 
by Lejeune et al. are based on the Kurucz (1995) library, extended to low temperatures 
based on empirical data and low-temperature stellar models: for M giants, spectra 
from Fluks et al. (1994) and Bessell et al. (1989,1991), and for M dwarfs, synthetic 
spectra from Allard & Hauschildt (1995), were used to correct spectral distributions with 
semi-empirical calibrations. There is also a version of the Y 2 set of isochrones using the 
color transformations by Green et al. (1987), but those transformed with the Lejeune et al. 
(1998) tables have been selected as those better reproducing the observed CMDs. 

The effects of the different color transformations are apparent in the [V, (V — R)} 
plane shown in the rigth panel. The Padova isochrones (Girardi et al. 2000) are slightly 
redder than the Y 2 ones along the unevolved main sequence and in the RGB, an effect that 
must be attributed to the color transformation alone since these parts of the isochrones 
match closely in the theoretical plane. The differences in the sub giant branch (SGB) 
region, mainly consisting in the Y 2 isochrones being brighter, are a combination of the 
differences in luminosity already noted in the theoretical plane, and the effects of the color 
transformation. In any case, the uncertainty in color is within the errors of the Hipparcos 
subdwarf data. 



6. ISOCHRONE FITTING TO THE CLUSTER CMDs 

Figure |7| displays a comparison among the fiducial lines of the three clusters. The 
cluster loci have been superimposed with no correction for possible differences in reddening 
or distance modulus. Note that, due to the presence of binary stars (see discussion in Paper 
II), drawing a fiducial sequence in the SGB region just above the main sequence turnoff is 
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not straightforward, especially for SL 556. Nevertheless, it is still evident from Figure [7| 
that the clusters show a trend in the V magnitude of the main sequence turnoff, and a 
(opposite) trend in the magnitude of the red-clump. The main sequence trend is clearly one 
in age, while the variation of the red-clump luminosity is consistent with the prediction by 
Girardi & Salaris (2001) about the luminosity evolution of the red-clump in the age interval 
1-3 Gyr. 

6.1. NGC 2173 

Figure |8] shows the comparison between the CMD of NGC 2173 and the Y 2 
(left) and Padova (right) theoretical isochrones. The same apparent distance modulus 
(m — M)v = 18.6 and metallicity Z=0.004, have been used in both cases, while due to the 
differences in the color transformations shown in Fig. || slightly different reddenings and 
ages are required for a good match with the two models. Both models reproduce quite well 
the main sequence, turnoff and SGB region, while they both fail to reproduce the color 
location of the RGB. A small change in slope would be needed to bring the theoretical Y 2 
RGB on top of the observed one, while a substantially different slope would be required to 
bring the Padova model into agreement with the data. In both cases a higher metallicity 
would improve the RGB fit, while the fit to the rest of the CMD remains similarly good. We 
adopt Z=0.004 for consistency with the two following papers in this series. The magnitude 
of the red-clump of the Padova models matches very well the observed one, while the color 
is slightly too blue, as all the AGB. Our crude estimate of the cluster age is therefore ~ 1.5 
Gyr. 

6.2. SL 556 

Figure |9| shows the comparison between the CMD of SL 556 and the theoretical 
isochrones from Y 2 (left) and Padova (right). A distance modulus of (m — M)y = 18.5 and 
a metallicity of Z=0.004 has been adopted in both cases, while slightly different reddening 
and age were used for the two models. The shape of the SGB region of this cluster is rather 
peculiar, and apparently quite different from theoretical predictions. As it will be shown 
in the next papers of this series, this shape can be explained by the presence of a fraction 
of binaries, also evident from the fact that there seems to be a second sequence, on the 
right, and parallel, to the cluster main sequence. Excluding this SGB region just above 
the turnoff, both models reproduce well the observed CMD. The Y 2 models also predict a 
correct RGB slope all the way up to the tip, while the Padova ones are steeper than the 
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observed RGB, for magnitudes brighter than the red-clump. Also in this case, the fit would 
be improved by adopting a higher metailicity, but then also a negative reddening would be 
required. 

6.3. NGC 2155 

Figure [It] shows the CMD of NGC 2155 compared with Y 2 (left) and Padova 
theoretical isochrones (right). Also for this cluster the same apparent distance modulus 
(m — M)v = 18.5 and metailicity Z=0.004, have been used in both cases, while slightly 
different reddenings and ages were required to achieve a good fit with the two models. Two 
isochrones are shown in the left panel, with or without core overshooting, the difference 
being only evident in the shape of the turnoff and SGB; some intermediate value of 
overshooting seems to be needed for a better fit of the observed points (see discussion in 
Paper II). The RGB is well reproduced by the Y 2 model all the way up to the tip, while the 
Padova one, due to the use of different transformations, starts to deviate from the observed 
points towards brighter magnitudes. A higher metailicity would improve the fit of the RGB 
slope, but then the other sequences would be too red in the model. The red-clump, only 
available in the Padova isochrones, perfectly matches the observed one, while the theoretical 
AGB is of course too blue, as the upper RGB. Using these models, our best estimate of the 
cluster age in this preliminary analysis is therefore between 2.5 and 3 Gyr. 

7. SUMMARY 

This preliminary investigation on three intermediate age LMC clusters, namely NGC 
2173, SL 556 and NGC 2155 is the starting point for a very detailed analysis aiming at 
a test on input physics for stellar models. In principle there is the opportunity to check 
the efficiency of convective core overshoot for star masses between 1.1 and 1.5 M & by 
comparison of these CMDs with synthetic CMDs and isochrones. The isochrone fitting 
to the CMDs of the clusters gives some information on their age and metailicity, but we 
must take into account that there are also uncertainties in the LMC distance modulus and 
interstellar reddening. Only a more refined analysis, considering also the star distribution 
in several regions of the CMD, will give reliable results. In the following two papers of 
this series, the synthetic CMD technique will be used independently by members of the 
Padova and Yale groups to determine the cluster characteristics, taking into account the 
uncertainties in the observations and in stellar evolutionary models, and to give some 
feed-back into the stellar evolutionary models themselves. 
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Fig. 1. SL 556 as imaged with F0RS1@VLT-UT1; North is up and East is to the left. 
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Fig. 3. — Statistical decontamination from LMC field stars of the CMD of NGC 2173. a) 
Contaminated cluster CMD, as measured in an annulus around the cluster center; b) Field 
CMD, as extracted from the outer part of the frame; c) Decontaminated cluster CMD; d) 
Stars subtracted from the CMD of panel a in order to obtain c. 




Fig. 4. — Statistical decontamination from LMC field stars of the CMD of SL 556. a) 
Contaminated cluster CMD, as measured in an annulus around the cluster center; b) Field 
CMD, as extracted from the outer part of the frame; c) Decontaminated cluster CMD; d) 
Stars subtracted from the CMD of panel a in order to obtain c. 
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Fig. 5. — Statistical decontamination from LMC field stars of the CMD of NGC 2155. a) 
Contaminated cluster CMD, as measured in an annulus around the cluster center; b) Field 
CMD, as extracted from the outer part of the frame; c) Decontaminated cluster CMD; d) 
Stars subtracted from the CMD of panel a in order to obtain c. 
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Fig. 6. — Comparison of the Y 2 (solid lines) and Padova isochrones (dashed lines) for 
Z=0.004 and ages 1, 2 and 3 Gyr. The theoretical plane is shown in the left panel, while the 
observacional [(V — R),V] CMD is shown on the right. 
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Fig. 7. — Comparison of the CMD loci for the three clusters. Note the sequence of ages in 
the main sequence turnoff. 
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Fig. 8. — Comparison between the CMD of NGC 2173 and theoretical models from the Yale 
group (left) and the Padova group (right). The adopted cluster parameters are listed in the 
figure labels. 
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Fig. 9. — Comparison between the CMD of SL 556 and theoretical models from the Yale 
group (left) and the Padova group (right). The adopted cluster parameters are listed in the 
figure labels. 
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Fig. 10. — Comparison between the CMD of NGC 2155 and theoretical models from the 
Yale group (left) and the Padova group (right). The adopted cluster parameters are listed 
in the figure labels. 
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Table 1. Log of observations 



Object Date RA DEC Filter Exptime Seeing 



NGC 2155 11 Jan 2000 05:58:32.4 -65:28:38.4 V 3768s -0.9 

R 2235s -0.9 

NGC 2173 12 Jan 2000 05:57:39.7 -72:58:41.4 V 928s -0.6 

R 855s -0.6 

SL 556 13 Dec 1999 05:32:25.0 -64:44:11.0 V 1856s -0.6 

R 1710s -0.6 
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